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It is common to relate the dynamic contact angle 0d to the relative speed between the 
substrate and the contact line; theory suggests 9^ oc U. In fact, available physical models 
show that the dynamic angle involves speed logarithmically and in a model dependent 
manner. Experimental data consistent with this interpretation is cited. 



1. Introduction 

One area of fluid mechanics that has been the subject of a large admixture of analysis, 
experiment and speculation is the subject of the moving contact line. A typical situation, 
common in many coating processes, refers to the contact line at the intersection of solid, 
liquid and gas regions, where the three-phase line moves relative to a solid substrate. A 
basic research question in this subject stems from the violation of the no-slip condition in 
the immediate neighborhood of the three-phase line of contact (e.g. Huh & Scriven 1971, 
for reviews see Dussan V. (1979), de Gennes (1985), and Kistler (1993)). As a result, 
within the usual continuum analysis, the stress diverges as the contact line is approached 
and the energy per unit length of the moving contact line is unbounded. This result may 
be viewed as an embarrassment of continuum modeling, but, in fact, it does indicate the 
need for a small cut-off length scale in macroscopic theories, as well as some more input 
from the physics at smaller length scales to properly interpret the meaning of any such 
cut-off scale. 

Perhaps the most basic feature of this problem is the aim to relate the local dynamic 
contact angle 9d(x), which is the arc tangent of the slope of the interface at a distance 
x from the contact line, to the local speed U with which the contact line moves over 
the substrate. For the case of a perfectly wetting fluid (vanishing equilibrium contact 
angle eq = 0) and small 8d, one finds 8^ oc U, which is known as Tanner's law. The- 
oretical justification for this result has been given (e.g. de Gennes 1985) and various 
generalizations have been offered. In dimensionless form, the speed is reported in terms 
of the capillary number C = Urj/^, which measures the relative importance of viscous to 
surface tension forces, where r\ is the fluid viscosity and 7 the interfacial tension. In fact, 
the functional form for the contact angle-speed relation is commonly written for small 
angles as Q\(x) w 9C ln(x / £ m i cro ) , where £ m icro is generally taken as a molecular length 
(e.g. Leger & Joanny 1992). Typically the capillary number varies over many orders of 
magnitude; values 10~ 7 < C < 10 _1 are common. The prefactor in this formula can be 
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important for interpreting experimental data, and so it is reasonable to interrogate more 
closely the functional dependence on speed. 

In this communication we wish to comment on one aspect of the moving contact line 
problem that has, perhaps surprisingly, been largely neglected and/or unappreciated. In 
particular, we note that detailed models for the perfectly wetting situation actually yield 
a dynamic contact angle versus speed relation 

9Ux)*9C\n(-^cA , (1.1) 



where (3 depends on the physical model introduced in the neighborhood of the contact 
line. We do not believe that it is necessarily appropriate to simply suppress the addi- 
tional dependence on speed (i.e. C) by replacing the argument of the logarithm by either 
tmacrol ^micro-, where these two lengths scales are taken as constants, or x/£ m i cro . Because 
of the large variation in C, not including this additional factor of capillary number when 



using (1.1) to interpret dynamical experiments may lead to significant discrepancies be- 



tween theory and experiment. Here we outline the basic idea behind (1.1) and present 
experimental evidence that supports the above interpretation. 

Another point that has received insufficient attention is the range of validity of equation 



(1.1). Near the contact line, (1.1) breaks down where x is of the same order as £ m icro- This 
restriction is evident as the general structure comes from a balance of viscous and surface 
tension forces alone. Not surprisingly, we estimate below that the microscopic scale is 
between several Angstroms and tens of Angstroms, depending on the microscopic forces 
assumed to be acting near the contact line. 

Towards large scales, x is commonly taken to be a static scale such as the capillary 
length or the size of a spreading drop (de Gennes 1985). Nevertheless, it should be noted 
that the flow near a moving contact line often resembles a coating flow, similar to the 
classical problem studied by Landau and Levich (e.g. Levich 1962). This leads to the 
appearance of another, dynamical length scale, that can become much smaller than the 
capillary length as the capillary number is small, which is typically the case. Thus a 



meaningful comparison between (1.1) and a macroscopic measurement of the dynamical 
contact angle might require a spatial resolution significantly below 1/10 or even 1/100 
of the capillary length. 

In the next section we will introduce two different models commonly used to treat 
moving contact line problems such as a spreading drop or a tape plunging into a pool of 
fluid. Then, in the third section we show that the lubrication equations corresponding 
to both models have similarity solutions for the interfacial shape that fix the functional 
dependence on the capillary number. The fourth section discusses the dynamical problem 
that equation ([O]) has to be matched to on an appropriate outer length scale. In the fifth 
section we explain measurable consequences of the two models for the dynamic contact 
angle, and discuss an experiment that helps to distinguish between them. We close with 
a summary and possible directions of future work. 



2. The model 

The usual dynamic balance for the steady flow "far" from the contact line involves 
capillary and viscous stresses. As the contact line is approached, the capillary-viscous flow 
leads to a stress singularity. A number of different physical effects have been suggested 
to relieve the singularity, and these either account for the fact that on very small length 
scales van der Waals forces act to maintain a finite thickness liquid layer on the solid 
substrate, or that at very high shear rates the boundary conditions and the transport 
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mechanism reference 

van der Waals Hervet & de Gennes (1984) 

Navier slip Huh & Scriven (1971) 

nonlinear slip Thompson & Troian (1997) 
shear thinning Gorodtsov (1990) 

diffuse interface Seppecher (1996) 

generalized Navier slip Shikmurzaev (1997) 

Table 1. Different models for the flow in the neighborhood of the contact line, with 

representative references. 



coefficients of the fluid are likely to be altered. Which model is appropriate might depend 
on the physical system at hand, or be a combination of the above. In Table [l] we provide a 
short overview of proposed physical models for flow in the neighborhood of a contact line; 
see also McKinley & Ovryn (1998). Most recently, there has been a considerable effort 
to base the understanding of the contact line physics on a microscopic, particle-based 
description, see for example Koplik et al (1989), Ruijter et al. (1999), and Abraham et 
al. (2002). The so-called "diffuse interface model" (see e.g. Seppecher (1996), Chen et 
al. (2000), and Pomeau (2002)) represents an intermediate approach, which models the 
liquid-gas interface as a Cahn-Hilliard fluid. This allows for example the extraction of 
effective interface equations (Pismen & Pomeau (2000)) different from those proposed by 
de Gennes's (1985). 

Below we restrict our attention to two different models which have proved particularly 
popular. The results are sufficient to highlight the measurable differences between differ- 
ent physical mechanisms. In model I, due to de Gennes and coworkers (e.g. Hervet & de 
Gennes 1982, de Gennes 1985), van der Waals forces are taken into account, so very close 
to the contact line there is a balance between surface tension and van der Waals stresses 
alone. In model II, proposed for example by Huh and Mason (1977) and Hocking (1977), 
the fluid is allowed to slip across the solid surface over a small slip length. 

For simplicity, we only consider the case of perfectly wetting fluids, i.e. of zero equilib- 
rium contact angle. Consistent with the local balances, the interface near the contact line 
remains nearly flat and we can use lubrication theory to describe the fluid motion. This 
approach amounts to a significant simplification of the mathematical treatment relative 
to the full two-dimensional flow problem (Cox 1986), but agrees with the full calculation 
when the dynamic contact angle is small. There are numerous indications that the small- 
angle theory in fact remains valid for slopes of order unity. For example, 0^/9 in equation 



(1.1) differs by only 2% from the full expression (Cox 1986), derived without the benefit 
of lubrication theory, up to a slope of 1. Also, Fermigier & Jenifer (1991) reported that 
small-angle theory holds experimentally up to an angle of 100°. 

To be able to describe an experiment like a flat plate plunging with velocity U into 
a reservoir of fluid (see Figure |l|), it is necessary to include other terms beyond the 
lubrication terms close to the contact line. Namely, we keep the full curvature term and 
include gravity. By doing this, the model is able to describe the crossover to a purely 
static, horizontal surface far away from the dynamical region. A basic unit of length is 
the capillary length £ c — which dictates the scale of the interface curvature far 

away from the contact line. In the van der Waals model I, which accounts for pressure 
variations owing to capillary, van der Waals, and gravitational forces, the equation for 
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Figure 1. A typical application involving a moving contact line: a plate plunges with velocity U 
into a liquid-filled container; the capillary length £ c — •v/-^- Since we assume wetting fluids, the 
meniscus creeps up the plate opposite the direction of motion and a nonzero dynamic contact 
angle 9d(x) is established. In the experiment of Marsh et al. (1993) referred to below the plate is 
replaced by a cylinder which can be tilted at different angles a. We also schematically indicate 
the different length scales relevant for this problem. 



(2.1) 



the stationary profile h(x) (cf Figure []]) is 

3C _ , 

where k is the curvature and a prime refers to differentiation with respect to x. Note that 
a positive C corresponds to the plate plunging into the fluid. A brief derivation of model 
I, as well as model II below, is given in appendix 1. The microscopic length parameter 
a, defined by 

a 2 = — 

67T7 

measures the strength of van der Waals forces relative to interfacial forces and is typically 
very small (on the order of Angstroms). 

Another distinct approach for the flow near the contact line is to introduce slip at 
the boundary, consistent with allowing the contact line to move parallel to the wall at a 
finite speed; the slip is a function of the shear rate. The simplest such law, introduced by 
Navier in the same paper that also enunciated the Navier-Stokes equation (Navier 1823), 
is 

du , 



(2.2) 



|j/=0 ~U — X—\y = 



(2.3) 



(see also Huh & Scriven 1971). Here U is the speed of the moving boundary, y = 
denotes the solid-liquid boundary, and A is the so-called slip length. A more complicated 
version of (2.3), in which A is itself a non-linear function of the shear rate, has been 
proposed in Thompson & Troian (1997). A standard calculation (appendix 1), leads to 
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the analogue of equation (2.1) for the slip model II 



3C 
h 2 



3A 



(1 



(2.4) 



The slip length A is usually considered to be on the order of tens of Angstroms. 



3. Scaling solutions near the contact line 

We now focus on the immediate neighborhood of the contact line, which we assume 
to be at x = 0. Owing to the flatness of the interface we can assume that n as —h" 
and gravitational influences can be neglect ed, b ut dynamical (viscous) effects have to be 
included. In the case of model I, equation ( |2.l| ) reduces to 

+ (3-D 



To make the dependence on parameters explicit, we note that (3.1) has the exact scaling 
solution 

h{x) = -£ 75 Ma- 1 C i/3 x), (3.2) 
where <f>x depends on the similarity variable £i = a~ 1 C 2 ^ 3 x and satisfies the equation 

*=_*f + 3§. (3.3) 

<l>\ 4>x 

Similarly, the lubrication approximation for model II gives 

S = -ti" ^ti" . (3.4) 
h z h 

In this case the scaling solutions are 

h(x) = A0 2 (A" 1 C 1/3 a;), (3.5) 
where the similarity variable is now £2 = \~ 1 C 1 / 3 x and the similarity equation is 



3 ±111 ^ ,„ 

<P2 02 



(3.6) 



Far away from the contact line in units of the microscopic lengths a and A, respectively, 
the solutions should be the same, resulting from a balance of classical viscous forces and 
surface tension. Indeed, as £ — ► 00 one finds to leading order 

&,a(0«3 3 / 3 £M#i,3)] 1/8 . (3.7) 

where the numerical constants b± t 2 have to be determined by numerical integration start- 
ing from the contact line. 

The boundary condition at the contact line incorporates the wetting behavior of the 
fluid. The basic assumption is that there is a microscopic length scale on which static 
forces dominate over dynamical ones, and a static profile can be assumed microscopically 
close to the contact line. 

In model I, following Hcrvet & de Gennes (1984), we are going to match (j>i(x) to 
a "maximal" film solution, corresponding to very strong wetting, whose thickness only 
goes to zero at (minus) infinity. However this maximal solution very closely approximates 
parabolic solutions of (|3.1|) that go to zero at some finite contact line position (Hervet & 
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de Gennes 1984). To leading order, we take the film solution to be of the form 

0i(a) = -i+eexp{e 1 3 /V3}, (3.8) 



and defer further details to appendix 2. Using (3.8) as an initial condition with adjustable 



parameter e, we integrate (3.3) towards £i — * oo. The parameter e is fixed to select the 
solution with van ishi ng curvature at infinity. Figure compares this sol utio n with the 
asymptotic form ( ^7 



We plot the rescaled slope tf>[ from the solution of ( |3.3| ) as the full 



curve and equation ([3/i]) with b\ = 1.44 as the dashed curve. This numerical value for b\ 
differs significantly from &i = 0.4 -3 1/6 « 0.48 given in de Gennes (1985), accounting for 
differences in normalization. We believe the difference is simply due to the large values 
of £ necessary for integration until a true asymptotic value is reached. 

Model II, on the other hand, can be extended down to 02(0) = 0, since the stress 
singularity was successfully removed. A vanishing equilibrium contact angle 9 eq — can 
thus be implemented by taking the boundary condition 02(0) = 02(0) = at the contact 
line and integrating (3.4) towards £2 - * 00. The corresponding value of 62 for the case of 
the Navier slip law was given in Hocking (1992). Thus the two constants, which establish 
the form of the interface profile, are 



h ~ 1.44, 6 2 = 3 1/3 exp(0.74/3) ^ 1.85. 



(3.9) 



4. Crossover to Landau-Levich-type behavior 



We now estimate the range of validity of the solution (3/7) as one moves farther away 
from the contact line. These ideas have close analogy to the classical analysis of Landau 
and Levich of a dynamical lubrication film (Levich 1962). Evidently, the contact line 
physics plays no role far away from the contact line, so the relevant lubrication equation 
is 

§ = -"'" <"> 

which is to be matched to a static profile at large distances. This equation has the general 
solution 

%)=4C ai+1/3 /(z/(4C ai )), (4.2) 
which has to be matched to a static meniscus on the capillary scale. This static solu- 
tion is characterized by a curvature that is approximately constant, thus h"(x) must be 
independent of C. This fact forces a\ = 1/3, so that we have 

h(x) = 4C 2 / 3 /(^/(4C 1/3 )), (4.3) 

which implies that the crossover will occur on a scale icC 1 ^ 3 , on which the logarithmic 
dependence in equation (3.7) for the slope begins to fail. 

The crossover to scaling of the form of equation (4.3) is demonstrated in Figure ||, by 
showing a full solution of equation (2.1), rescaled according to (3.2). Results are given 
for two (small) values of the capillary number differing by a factor of ten. Again, the 
free parameter in the maximal film solution d3.8| ) is used to shoot for the flat interface 
corresponding to the surface of the fluid-filled container. For small values of £1, the 
solution corresponds to the lubrication form given before, while on a scale x/£ c f=a C 1 / 3 
the transition to the Landau-Levich region is observed. In rescaled coordinates £1 the 
location of this crossover should thus be proportional to C itself, as is clearly seen from 
Figure 2. We have chos en t he smaller of the two values of C such that the region over which 
the asymptotic form ( |3.7| ) of the interface can be applied is zero, to highlight possible 
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problems in comparing asymptotic solutions with experimental data. To interpret the 
measured dynamical contact angle equation (3.7) is no longer sufficient, but the full 
solution of the similarity equation (3.3) has to be considered. 



5. Comparison with experiment 

It is common practice in the literature (e.g. de Gennes 1985, Cox 1986) to consider the 
derivative of the profile h(x), evaluate it at some macroscopic distance from the contact 
line x = £ m acro, and to interpret the slope of the interface in terms of the so-called 
"dynamical contact angle", taxi9d(x) — dh/dx. T his a pproach is the common o ne t aken 
i n ex periments as well. Thus, using the solution (3/7) in the similarity forms (3.2) and 
( |3.5| ) for models I and II, respectively, and neglecting lower-order terms, we obtain 

6 3 dyn (x) =9Cln(£ mocro /L 1 , 2 ), (5.1) 

where are microscopic lengths appropriate for each model. There are two fundamen- 
tal issues with this approach: First, depending on the experimental system it is not clear 
what is the best choice for £ ma cro- Second, what is usually taken as a fixed microscopic 
length Li,2 is actually strongly dependent on the capillary number. Namely, the two 
models give 

L x = aC-^/h, and L 2 = ACT 1/3 /V (5.2) 

as established in (3.5) and (3.5), respectively. 

In particular, the C dependence that appears in the microscopic length is different in 
the two models. It is also clear that it is impossible to interpret L± t 2 directly in terms 
of some fixed microscopic length near the contact line, but rather it is a dynamical 
quantity. To our knowledge, this fact has never been appreciated in either theoretical or 
experimental work. This observation appears to be significant, since by comparing the C- 
dependence it potentially allows one to distinguish between different microscopic models 
from a macroscopic measurement. Chen and Wada (1989) imaged the profile near the 
con tact line of a spreading droplet and so provided the first experimental confirmation 
of (|]l]). However, owing to the small range of capillary num bers studied, it is difficult 
to distinguish between the two lengths defined in (|5.2| ). Below we will therefore 
concentrate on another experiment (Marsh, Garoff & Dussan 1993), which allowed C to 
be varied over more than two orders of magnitude. 

It is also interesting to note that the logarithmic dependence on capillary number was 
only obtained in the fully nonlinear treatment outlined above. In the classical studies of 
the flow in the neighborhood of the dynamic contact angle, for example in Cox (1986), 
a matched asymptotic analysis is used which further assumes a form involving integer 
powers of C and does not recognize that the scale of the inner region can itself involve 
C. However the form of the asymptotics ( |5.l|) , implying a logarithmic dependence on C 
in the full solution, clearly shows that the profile cannot be expanded in integer powers 
of C. As a result, the classical analyses are not able to identify the sort of dependencies 
given by ( |5.2| ). Although these dependencies are only logarithmic, as mentioned above, 
the capillary number often varies through many orders of magnitude in experiments, so 
the logarithm in (|]l]) can in general not be approximated by a constant, as is most often 
done (e.g. King 2001). 

In a recent experiment, the effect of large variations of C on the contact line was 
investigated very carefully by Marsh et al. (1993), who measured the dynamic contact 
angles on a cylinder plunging at an angle into a liquid bath. These authors essentially 
used the form (5.1) to fit the whole shape of the interface close to the contact line, and 




Figure 2. The rescaled slope <fi[ = h'a^ 1 C 1/:i versus log 10 (£i) = \og 10 (a~ 1 C 2 ^ 3 x 
C = 10" 5 ,10" 6 in the geometry of Figure |l[ The ratio a/l c equals 10" 8 . The full line _ 
solution of (3.1), the dashed linecorresponds to the asymptotic form (|3.7|). The two dotted Unes 
are solutions of the full system ( p.l[ ) including gravity, marked with their values of the capillary 
number. Note that for the smaller capillary number the size of the overlap region where the 
asymptotic form (3.7) can be applied has shrunk to zero. 



foi- 
ls a 
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Figure 3. A plot of the microscopic length L, taken from Figure 4 of M arsh et al. (1993) 



function of the speed U. Their L is the equivalent of Z/1,2 as given in (5.1). The angle a refers to 
different tilts of the solid relati ve t o the liquid surface. To make a comparison with the dynamical 
length sales Li,2 as given in j5.2| ), we have added to the figure the solid and the dashed lines 
with slope 2/3 and 1/3, respectively. 

included static contributions to account for the effects of surface tension and gravity 
away from the contact line. (Note that instead of the third power on the left hand side, 
they actually used a more complicated function g(x), but which becomes g(x) ps x 3 /9 
for small arguments. This limit is relevant for the small angle case we are studying here.) 
This approach leaves out dynamical effects of the kind predicted by Landau and Levich 
(Levich 1962), which are important on an intermediate scale between the microscopic 
ones and the capillary length, and should be taken into account in a more refined theory. 
Marsh et al. (1993) treat the static contact angle 9 eq (called 9 act by the authors) as a free, 
and possibly C-dependent, parameter, to be determined from experiment. The authors 
conclude that for their system 9 eq = 0, which is the case treated here. 

From the fit of (|5.l|) to their data, Marsh et al. (1993) extract a length L, which 



is found to depend significantly on capillary number, as suggested by (p.2[)! They also 
report L to be independent of the tilt angle a within experimental error, which further 
emphasizes that the response is dominated by local features. In Figure || we present a 
plot of the measured length as function of capillary number, and compare it with the 
slopes suggested by the van der Waals model I and the Navier slip model II, respectively. 
Although it is difficult to draw firm conclusions owing to the large scatter in the data, 
2/3 seems to be favored. Using the two different fits plotted in Figure 0, we are also 
able to determine the cut-off lengths a and A, assuming that the corresponding physical 
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mechanism is really relevant for the particular materials involved. We find a ~ 4Afor the 
van der Waals model and A ~ 30Afor the slip model. Using the value of A = 10 _20 J 
(Russel, Saville & Schowalter 1989) for the Hamaker constant for water and an adjacent 
solid surface, and 7 = 0.07N/m, we find a ~ lA, consistent with the above value. 
However, it is important to keep in mind that there is no reason why a single cutoff 
mechanism should necessarily dominate in the experiment, which would lead to still other 
exponents. Additional mechanisms for relieving the contact line singularity are listed in 
Table ^ which is the dominant mechanism could also depend on capillary number and, 
in particular, on the type of solid substrate or fluid involved. 



6. Conclusions 

We have considered flow local to a moving contact line using a lubrication approach. 
Our basic message has been to indicate that the appearance of logarithmic corrections 
in capillary number to the usual "Tanner's law", 8^ <x C, are a general feature of the 
mechanical response. The interpretation of the results is that the "microscopic" length 
scale that is involved when supplying a small-scale cut-off to relieve the well-known 
stress singularity in the moving contact line problem is a dynamical (speed dependent) 
quantity. Experimental data consistent with this interpretation is cited, and microscopic 
parameters that come from the comparison with different theories agree with physical 
considerations. 

Since important applications of contact line theories apply to angles up to 180°, it 
would be very useful to extend the lubrication theory considered here to a full two- 
dimensional treatment of the flow in the corner region. We suspect that large parts of 
the calculation in Cox (1986), which erroneously assumes a regular expansion in C, could 
be adapted to a proper similarity description. This means that the dynamical contact 
angle has the scaling form 9d = f(C, x/Lx^), where L\ i is one of the dynamical length 
scales defined by (5.2), with corresponding forms for the velocity field as well. Such a 
theory might be able to explain more recent experiments (Chen, Rame & Garoff 1995) on 
moving contact lines performed at higher capillary numbers, yielding dynamical contact 
angles of up to 155°. 

Furthermore, the velocity of the contact line relative to the substrate is in general 
not 'perpendicular to the contact line, as highlighted in recent experiments of droplets 
running down an inclined plane (Podgorski, Flesselles & Limat 2001). In this case the 
flow is truly three-dimensional, and it may no longer be sufficient to simply project the 
velocity onto the normal to the contact line (Blake & Ruschak 1979). Such a three- 
dimensional description would be necessary to complete our understanding of corner 
singularities that form at the back of running drops (Stone et al. 2002), and may apply 
to a range of other contact line phenomena as well. 



7. Appendix 1: Derivation of lubrication equations 

Here we briefly recall the derivation of the interface, or lubrication, equations for thin 
viscous films (Levich (1962)). For pressure-driven flow along the surface and absorbing 
hydrostatic pressure variations into the pressure p, the velocity parallel to the plate can 
be represented as a second-order polynomial 

u = a + a x y + y 2 ^- , (7.1) 
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where y is the distance normal to the plate. At the free surface y = h(x), shear gradients 
du/dy vanish, giving a± = —p'h/rj. Finally, from the slip condition ( [2.3] ) we have ao = 
U - Xp'h/r). 

Since the contact line is stationary, the mass flux through the film is zero everywhere, 
/ u(y) dy = 0, and thus 

= (U - \p'h/r))h - p'h z /iri (7.2) 

is the equation for the film profile h(x). In the presence of van der Waals forces, the 
dynamic pressure in the liquid is 

A 



p = 7« 



6irh 3 



gpx, 



(7.3) 



where A is Hamaker's constant. Substituting (7.3) into (7.2) and assuming A = gives 
(O), while A = at finite A leads to (O). 



8. Appendix 2: expansion for the maximal film 

Here we give some more details on the solution of (3.1) for the "maximal film" of 
Hervet and de Gennes (1984). The general form of the film profile is 



-. oo 

ba - 7 2^ 767' 

s i=0 S 



(8.1) 



where wc denote the similarity variable by £. This expansion has no free parameters, as 



the values of the coefficients a,; are obtained directly from substituting (8.1) into ( J3 . l[ ) . 
We find 



a 



T , a\ 



-2/5, a 2 = -1764/275, . 



•2) 



However, there is a one-parameter family of solutions of (3.1) tha t decay for £ — > — oo. 
This solution is found by linearizing around the base solution (8.1), i.e. = <f>o(£) + 

my- 

(5(60o + 4<t>o<P'o) - 35' + S"% = 0. (8.3) 
Equation ( |S.3| ) is solved using a WKB-type ansatz, 

e 



5{£) = eexp 



V3 



(8.4) 



The 0(£°) contribution in the exponent turns out to be a logarithm, so the full struc 
ture is 

' " he- 3 *] 



72 ex P i J2 



V3 



and the coefficients are found to be 

b = 1 , b 2 = , 6 3 = 32/15 , 64 = 9V3/5, . 
Thus the general form of the solution in the film region is 



5) 



(8.6) 



■>) 



with a single free parameter e. An alternative description would be an expansion of the 
form 



*(0 = E§. 
1=1 15 
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with c\ = — 1 and C2 a free parameter. However, the convergence of the asymptotic series 



(3.8) turns out to be very bad, as perhaps is to be expected from the structure of the 



WKB solution. 
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